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Abstract 

We summarize results of recent studies of heavy quarkonia correlators and spectral functions at finite temperatures 
from lattice QCD and systematic T-matrix studies using QCD motivated finite-temperature potentials. We argue that 
heavy quarkonia dissociation shall occur in the temperature range 1.2 < T c ilT c < 1.5 by the interplay of both screen- 
ing and absorption in the strongly correlated plasma medium. We discuss these effects on the quantum mechanical 
evolution of quarkonia states within a time-dependent harmonic oscillator model with complex oscillator strength and 
compare the results with data for 7?aa/^aa M fr° m RHIC and SPS experiments. We speculate whether the suppres- 
sion pattern of the rather precise NA60 data from In-In collisions may be related to the recently discovered X(3872) 
state. Theoretical support for this hypothesis comes from the cluster expansion of the plasma Hamiltonian for heavy 
quarkonia in a strongly correlated medium. 
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1. Introduction 

In Physics we know fields of research which are 
opened by ingenious experimental discoveries and 
sometimes have to wait even decades for their proper 
theoretical explanation. This was so in the case of su- 
perconductivity, where after Kamerlingh-Onnes' obser- 
vation of the vanishing resistance in Mercury at liquid 
He temperatures |l]] more than four decades of theory 
development were necessary before a satisfactory ex- 
planation could be formulated by Bardeen, Cooper and 
Schrieffer 0), Sometimes it is vice-versa, like in the 
case of Bose-Einstein condensation [3]. The verifica- 
tion of the theory formulated by Bose and Einstein in 
1924/25 succeeded only 70 years later J3l since ex- 
tremely subtle experimental techniques had to be devel- 
oped to cool a sufficient number of atoms in a trap to 
nano-Kelvin temperatures. 

The effect of J/if/ suppression (more general, heavy 
quarkonium suppression) was suggested by theory 15J 
to be a signal of quark-gluon plasma formation before it 
was actually seen in first heavy-ion collisions at CERN 
SPS and thus at first glance seems to belong to the lat- 
ter class of discoveries. However, after 24 years of in- 
tense experimental research and theory developments, 
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one has the impression that heavy quarkonia became a 
field where theory and experiment have to work hand 
in hand in order to make progress, like at this work- 
shop. A key problem to be solved is that the usage of 
heavy quarkonia states as a probe for the diagnostics of 
the quark-gluon plasma (QGP) state of matter in ultra- 
relativistic heavy-ion collisions requires the knowledge 
of a baseline, e.g., from their production and evolution 
characteristics in situations when a QGP is absent. For 
a recent review see, e.g. |0] and references therein. Cur- 
rent issues in the experiment-theory dialogue are sum- 
marized, e.g., in 10]. Aspects of quarkonium production 
at LHC are discussed, e.g., in [8, 9]. For these proceed- 
ings we first summarize basic theory issues on the use 
of quarkonia as probes of the QGP before outlining the 
aspect of the modification of quarkonium formation by 
a strongly correlated QGP within a quantum mechanical 
model. 

The spectrum of low-lying heavy quarkonia states is 
perfectly decribed by solutions of the Schrodinger equa- 
tion for confining potentials of the Cornell-type II 1011 . 
While at zero temperature the Cornell potential is nicely 
reproduced by lattice QCD simulations for the change in 
the singlet free energy of the system when static color 
charges are inserted, it has been argued that in a ther- 
mal system the internal energy should be used as the po- 
tential instead 111]. The proper form of potential to be 
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used for studying quarkonia states in the medium may 
be a superposition of both forms depending on the in- 
terplay of time scales for interaction and thermal relax- 
ation in the heat bath [12]. The status of this discussion 
is unsettled. In the literature (cf. II 1311 ) calculations with 
both potential models derived from fits to lattice data 
are found whereby the use of the free energy as a heavy- 
quark potential implies a weaker binding and thus lower 
temperatures for the Mott transition of heavy quarkonia 
states, see also 114 1. 



come the threshold for impact dissociation of quarkonia 
1 17, 18 , 1^1 . see the lower panels of Fig. Q]for corre- 
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Figure 1: The "classical" picture: modification of binding energies 
of heavy quarkonia states in a hot plasma by static screening leads to 
their dissociation at the corresponding Mott-temperatures. 



In Fig.Q]we illustrate the "classical" picture of bound 
state dissociation in a hot plasma by screening of the in- 
teraction for the example of a screened Cornell potential 
applied to charmonium and bottomonium states IU5I1 : at 
state dependent (Mott) temperatures, the binding energy 
vanishes and bound state merges the continuum of scat- 
tering states. Two comments of significance for experi- 
ments are in order: (i) just above the Mott temperatures, 
remnants of the charmonia states survive as resonances 
in the continuum to be identified by peaks in the spec- 
tral functions and strong correlations in the correspond- 
ing scattering phase shifts 11611 : (ii) before reaching the 
Mott temperatures, the binding energies are already suf- 
ficiently lowered so that collisions with particles from 
the medium may have sufficient thermal energy to over- 



sponding estimates. This effect is particularly dramatic 
for T which may undergo thermal dissociation already 
at T ~ 250 MeV while its Mott temperature is ~ 470 
MeV. 

Both effects tend to wash out the pattern of sequential 
suppression for heavy quarkonia states expected from 
the "classical" picture of the Mott effect ||20j]. The de- 
scription of quarkonia states in a hot QGP medium in 
the vicinity of the critical temperature should therefore 
treat bound and scattering states on an equal footing. 
This is appropriately achieved within a thermodynami- 
cal T-matrix approach, which has been developed to ad- 
dress the spectral properties of quarkonia 112 ill as well 
as open flavor meson states II221I23I1 . Such an approach 
allows a simultaneous description of quarkonia suppres- 
sion and heavy flavor diffusion, presently under scrutiny 
for RHIC experiments, see 1 1311 and references therein. 

The lowering of the quarkonia binding energy by 
screening has important consequences also for the pro- 
cess of quarkonium formation itself, once it occurs in a 
QGP environment. We will illustrate this aspect within 
a simple model. 

2. Time-dependent harmonic oscillator model 

For our discussion of the quantum mechanical evolu- 
tion of quarkonia in an evolving QCD plasma state, we 
will employ here a generalization of the harmonic os- 
cillator model 12411 to time-dependent one with complex 
oscillator strengths (THO model). Aspects of an opti- 
cal potential for the propagation of charmonia through a 
medium have been already discussed, e.g., for cold nu- 
clear matter in 1 25 , 26| and for a quark-gluon plasma 
in 11271 l28ll . The merit of such a model is its simplicity 
and transparence as well as tractability. We consider the 
time-dependent Hamiltonian for heavy quarkonia in the 
form 



H(t) = 2m Q + ?- + W(T)r 2 (r) , 

2yU 2 



(1) 



where fi = m^/2 is the reduced mass and niQ the 
heavy quark mass. The complex oscillator strength 
o> 2 (t) = u) 2 r (t) + icj 2 (t) has an implicit time dependence 
due to the temperature evolution T(r) of the system sur- 
rounding the evolving heavy quarkonium state. The 
temperature dependence of cl>r(T) resembles screening 
or strengthening of the quarkonium interaction, while 
a nonvanishing a>i{T) signals for quarkonium absorp- 
tion or dissociation processes. The quadratic depen- 
dence of the imaginary part of the (optical) oscillator 
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potential is motivated by the phenomenon of color trans- 
parency, see also l29ll and references therein. The con- 
ditions w/(0) = and wr(0) = &ty apply for the vacuum, 
where the spectrum of the low-lying quarkonium states 
is approximated by a suitably chosen oscillator strength 
co = City = const. 

The general classical trajectories for the Hamiltonian 
(Q} are linear combinations of the two solutions 



r(f)=p(f)exp(±/0(f)), 0( f ) 



=1 



df 



The amplitude p(t) fulfills the Ermakov equation 



p(t) + oj z (t) p(t) 



1 



p\t) 







(2) 



(3) 



for which exact solutions exist 1 30], allowing to evaluate 
the time evolution operator using path integral methods 

1/3 



U(tf,ti) = 



VP/Pi 



2ni sin(0y - </>,) 



(4) 



where the classical action functional S c \[p{t)\ enters. 
The survival probability (suppression factor) for i/ifi as 
defined in [24] can be generalized to the THO case and 
applied to the QGP diagnostics in collisions of heavy 
ions with mass number A when identified with the ex- 
perimentally determined quantity 



Raa 

^CNM 



Pf/Pi 



(5) 



where accounts for the cold nuclear mat- 

ter (CNM) effects from charmonium absorption 
and modification of charm production by shadow- 
ing/antishadowing of gluon distribution functions in the 
CNM of the colliding nuclei. Both effects have been 
discussed in contributions to this workshop and are in 
principle accessible by analysis of pA collision exper- 
iments, see l32ll and references therein. We restrict 
our discussion here to ground state charmonium at rest 
in the QGP medium (pj — y = 0) so that the discussion 
of Lorentz boost effects on the formation process can 
be omitted and also a detailed discussion of feed-down 
from higher charmonia states will be given elsewhere 

In the following we show that the anomalous J /iff 
suppression in both, SPS and RHIC experiments can 
be simultaneously described with the natural assump- 
tion that above the citical temperature relevant screen- 
ing parametrized with a temperature dependent, com- 
plex oscillator strength. The time evolution of temper- 
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Figure 2: Recent results (data points) for the equation of state from 
lattice QCD j34j] compared to the fit formula employed here. 



PHENIXAuAu,T =3ll MeV 
PHENIX AuAu,T a = 199 MeV 
NA6I1 PbPh, T fl = 224 MeV 
NA60 PbPh, T = 190 MeV 




6 8 10 
x [fm/c] 



Figure 3: Temperature evolution across the QCD transition for differ- 
ent initial conditions in SPS and RHIC experiments. 



ature itself will be given by longitudinal (Bjorken scal- 
ing) hydrodynamic evolution of a fireball volume V(t) 
under entropy conservation 

S(t) = const = s(T(i)W(t); V(t) = A T z(t) (6) 

with initial conditions determined by a Glauber model 
for nucleus-nucleus collisions. The temperature de- 
pendence of the entropy density s(T) is taken from 
recent lattice QCD simulations 13411 which are well 
parametrized by the simple ansatz 



s(T) = 9.0 T 3 



1 + tanh 



( T-T \ 
1 0.54 7V 



(7) 



with 7q = 0.189 GeV, see Fig. [2] The resulting temper- 
ature evolution is shown in Fig. [3] for initial values of 
temperature which correspond to initial entropy densi- 
ties given by the Bjorken formula l35ll so = ^7" L=o' 
where the transverse overlap area At is taken from 113611 . 
For the initial time to of the thermodynamical fireball 
evolution is assumed to depend on the center of mass 
energy of the collision and we take here 0.6 (1.0) fm/c 
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for RHIC (SPS) experiments. This completes the def- 
inition of the the THO model for applications to QGP 
studies in heavy-ion collisions which we discuss in the 
next section. 

3. Anomalous suppression and the In-In case 

A key indicator for QGP formation in heavy-ion col- 
lisons is anomalous suppression, the deviation of exper- 
imental data for the J/ip production ratio (0 from unity. 
This effect, first observed at CERN SPS for Pb-Pb col- 
lisions at = 17 GeV, has been qualitatively con- 
firmed by RHIC experiments with Au-Au interactions at 
V? = 200 GeV whereby surprising new findings were 
revealed: (i) the suppression is stronger at forward and 
backward rapidities rather than at midrapidity where the 
particle densities are highest, (ii) the onset of anomalous 
suppression and its dependence on centrality scales with 
the charged particle density at midrapidity rather than 
with energy density. While (i) is caused mainly by an- 
tishadowing and to some extent by geometry 13711 . the 
second finding is still not understood. A third puzzling 
issue raised by Carlos Lourencp in discussions at this 
meeting is (iii) the dip in the centrality dependence of 
the J/if/ suppression ratio which is seen in the rather pre- 
cise data of the NA60 collaboration for In-In collisions 
and so far widely ignored by theorists. 

We want to report rather fresh results within the THO 
model which might get substantially improved in near 
future but the main idea of which is in the spirit of this 
workshop, namely to provide a more appropriate the- 
oretical basis for future experiments on heavy quarko- 
nium production, in particular in the LHC era. We 
are convinced that any deeper insights into the quan- 
tum mechanical evolution of the cc state in the hot 
QGP medium can become very important in this con- 
text. Here we will apply the THO model to extract the 
real (imaginary) part of the oscillator strength and thus 
the amount of screening (absorption) as a function of the 
temperature of the QGP medium from the experimental 
data. 

In Fig. |4] we present a fit which ignores the In-In dip 
and is governed by the screening of the confining inter- 
action which monotonously drops to zero in the temper- 
ature range T/T c ~ 1.2 - 1.6, in gross accordance with 
the analysis of the charmonium spectrum from the tem- 
perature dependent potential on the basis of the heavy- 
quark free energy from lattice QCD simulations with a 
Mott temperature below 1.4 T c . This fit is rather insen- 
sitive to strong variations of the absorptive part of the 
potential. For recent systematic, quantum field theoretic 
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Figure 4: Anomalous suppression within the THO model com- 
pared to data from NA50, NA60 and PHENIX (upper panel) and the 
temperature dependence of the real and imaginary parts of the HO 
frequency (lower panel). 



motivation of complex charmonium potentials see, e.g.. 
Refs. 



In Fig. [5] we present a possible parametrization of the 
complex THO model accounting for the dip in the In- 
In data liJJ, 41] while still being compatible with the 
NA50 Pb-Pb and the PHENIX Au-Au data which have 
considerably larger error bars. The dip reflects a non- 
monotonous temperature behaviour of the confining po- 
tential due to a resonance-like contribution to the oscil- 
lator strength which results in a J/i/r regeneration pat- 
tern as a function of the charged particle multiplicity. 
This parametrization is more susceptible to changes in 
the absorptive part in the complex oscillator strength: if 
a>i(T) ~ 0)^/2 a too strong J/t/f suppression would re- 
sult. 

Let us present a speculation about the possible origin 
for a resonant-like strengthening of the heavy-quark po- 
tential in the vicinity of the chiral/ deconfinement transi- 
tion. The Lippmann-Schwinger equation for p-J/tf/ scat- 
tering in the D-D* channel can be depicted in the fol- 
lowing way 
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Figure 5: Same as Fig. [4] but with a temperature dependence of the 
HO frequency which exhibits a resonance-like structure (lower panel), 
allowing for a fit of the "dip" in the In-In data (upper panel). For 
details, see text. 
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and its solution may produce a new state or resonance, 
such as the X(3872) recently discovered at Belle and 
BaBar. For a recent discussion, see 0.52], and refer- 
ences therein. The scattering kernel, 
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is given by double quark exchange. It couples both sides 
of the X(3872) medal: the D-D* side which has led 
to the hypothesis it may be a molecule made of these 
states, and the p-J/iff side, the channel to which it ac- 
tually predominantly decays. Here we want to make 
the contact with plasma physics where a similar process 
contributes to the plasma Hamiltonian for two-particle 
states in a strongly correlated, nonideal plasma. It is 
represented diagrammatically as 
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and follows from a cluster expansion of the plasma 
medium, for details see Ref . ll43ll . In the present context, 
we may identify the T-matrices with resonant mesonic 
states: T = p, T = D and T" = D*. The resulting 
contribution to the plasma Hamiltonian for quarkonia 
states is proportional to the partial density of p mesons 
in the medium which is peaked just above the QCD 
transition temperature, so that a contribution with the 
shape of <jl>i,(T) in the lower panel of Fig. [5]could be ex- 
pected. A more detailed calculation could be based in 
Refs. 1 44, 45] and is in progress J33ll . Quark exchange 
processes for quarkonium suppression have first been 
introduced within nonrelativistic potential models 145 ] 
and successfuly been applied to describe CERN SPS 
data 14611 . Their reformulation within relativistic quark 
models has confirmed results for the the behaviour and 



magnitide of the J/if/ dissociation cross section by n (47 ] 
and p meson impact 14811 . These models can be used 
to derive formfactors for the otherwise very powerful 
chiral Lagrangian approaches to charmonium dissocia- 
tion, see 14911 and references therein. In the context of 
the X(3872) conjecture the finding of Ref. ifsoh is im- 
portant that the spectral broadening of D-mesons due to 
their Mott effect at the chiral transition entails a quali- 
tative increase in the dissociation rate by meson impact. 
Here the p meson plays the dominant role lf5ltl . which 
motivates the conjecture that the ladder-type iteration of 
the quark exchange interaction kernel Um P may provide 
sufficient strength to produce the X(3872) state in that 
channel. 

4. Conclusion 



Recent studies of heavy quarkonia correlators and 
spectral functions at finite temperatures in lattice QCD 
and systematic T-matrix approaches using QCD moti- 
vated finite-temperature potentials support that heavy 
quarkonia dissociation shall occur in the temperature 
range 1.2 < T ( i/T c < 1.5 whereby the interplay of 
both screening and absorption processes is important. 
We have discussed these effects on the quantum me- 
chanical evolution of quarkonia states within a time- 
dependent harmonic oscillator model with complex os- 
cillator strength and compared the results with data for 
Raa/Raa(CNM) from RHIC and SPS experiments. Be- 
sides the traditional interpretation, with a threshold for 
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the onset of anomalous suppression by screening and 
dissociation kinetics at dN c h/dj] ~ 300, we suggest 
an alternative arising from the attempt to model the 
dip the suppression pattern of the rather precise NA60 
data from In-In collisions at dN c h/dr] ~ 150 - 250. 
We suggest that this dip indicates the true threshold 
for the onset of anomalous suppression due to the cou- 
pling of charmonium to the DIDq channel with the re- 
cently discovered X(3872) state. Although details need 
to be worked out, the theoretical basis for supporting 
this hypothesis has apparently been developed in plasma 
physics with the concept of a plasma Hamiltonian for 
nonrelativistic bound states like heavy quarkonia when 
they are immersed in a medium dominated by strong 
correlations like bound states, to be systematically ad- 
dressed within cluster expansion techniques. This illus- 
trates that in the studies of the interrelation of quarko- 
nia with the QGP there are still many interesting and 
challenging aspects to be clarified within the further de- 
velopment and close collaboration of theory and exper- 
iment. 
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